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Abstract 

l_H ■ We show that given an ordinary differential equation of order four, it may be possible to 

determine a Lagrangian if the third derivative is absent (or eliminated) from the equation. 
This represents a subcase of Fels'conditions [M. E. Fels, The inverse problem of the calculus 

* *""] ' of variations for scalar fourth-order ordinary differential equations, Trans. Amer. Math. 

Soc. 348 (1996) 5007-5029] which ensure the existence and uniqueness of the Lagrangian 
in the case of a fourth-order equation. The key is the Jacobi last multiplier as in the case 

_ , of a second-order equation. Two equations from a Number Theory paper by Hall, one of 

^ ' second and one of fourth order, will be used to exemplify the method. The known link 

between Jacobi last multiplier and Lie symmetries is also exploited. Finally the Lagrangian 
of two fourth-order equations drawn from Physics are determined with the same method. 

^ ; 

1 Introduction 

oo ! 

It is well-known that a Lagrangian always exists for any second-order ordinary differential 
equation [30]. What seems less known is that the key is the Jacobi last multiplier [TT], |12| . 
|13j . [30] . which has many interesting properties, a list of which can be found in |26j. 

It is a matter for historians alike to find out why Darboux [3], Helmholtz [10], Koenigs- 
berger himself] and many other successive authors, e.g. Douglas [I] and Havao [7J, never 
acknowledged the use of the Jacobi last multiplier in order to find Lagrangians of a second- 
order equation. 



'Permanent address: Department of Mathematics, University of York, Heslington, York YO10 5DD, U.K., 
e-mail: ama5@york.ac.uk 

x In 1902-1903 Koenigsberger wrote Helmholtz's biography [T7] - which in 1906 was (abridged) translated into 
English with a Preface by Lord Kelvin [19] - after he wrote his 1901-book on Mechanics [16]. Neither books cite 
the connection between Jacobi last multiplier and Lagrangians. In 1904 Koenigsberger wrote Jacobi's biography 
[18] where the Jacobi Last Multiplier is extensively described. 

2 Havas even cites the book by Whittaker |30] but only in connection with the formulation of Lagrangian 
equations. 
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The Norwegian Sophus Lie, who carefully studied Jacobi's work [9], found a connection 
between his groups of transformations and the Jacobi last multiplier [20], [21]. The Italian 
Bianchi presented Lie's work in his lectures on finite continuous groups of transformations, 
and described quite clearly the Jacobi Last multiplier and its properties pp. Neither of them 
cited the connection with Lagrangians. 

A deluge of papers have been dedicated to the solution of the inverse problem of calculus 
of variation, namely finding a Lagrangian of differential equations. Most of them did not 
acknowledge the seminal work of Jacobi. Some did (see references in [26], and [28]). 

In [5] Fels derived the necessary and sufficient conditions under which a fourth-order equa- 
tion, i.e. 

= F{t,u,u',u",u"'). (1) 

admits a unique Lagrangian, namely: 

d 3 F 
d{u'"Y 

dF 1 d 2 / dF \ d / dF \ 3 dF d / dF \ 1 dF dF 1 / dF \ 3 
cV + 2dfi ~ dt \dv") ~ 4 9u^dt [dv 7 " J + YdvP^F' + 8 \dvF' ) 

Here we propose to extend the use of the Jacobi last multiplier in order to find the Lagrangian 
for ordinary differential equations of order four satisfying a subset of Fels' conditions ((SI)-©. 

The paper is organized in the following way. In the next section we present the properties of 
the Jacobi last multiplier and its connection to Lie symmetries; then we show the connection 
of the Jacobi last multiplier with Lagrangians for any second-order equation, and how and 
when can this connection be extended to fourth-order equations. In section 3 two equations 
from a Number Theory paper by Hall [6j, one of second and one of fourth order, are used to 
exemplify the method. The link between Jacobi last multiplier and Lie symmetries [20], |21j 
is also exploited. In section 4 we consider two fourth-order equations drawn from Physics and 
use the method of the Jacobi last multiplier to find their respective Lagrangians, and lastly 
we provide some of the many examples where the method does not work. The last section 
contains some final remarks. 



2 Jacobi Last Multiplier 

The method of the Jacobi Last Multiplier [H]-[13j) provides a means to determine all the 
solutions of the partial differential equation 

A/ = f>(xi,...,*n)J^=0 (4) 
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or its equivalent associated Lagrange's system 

dxi dx2 
ai a 2 



dx n 



(5) 



In fact, if one knows the Jacobi Last Multiplier and all but one of the solutions, then the last 
solution can be obtained by a quadrature. The Jacobi Last Multiplier M is given by 



9(f, U\,U 2 , ■ ■ ■ ,^n-l) 

d(x 1 ,x 2 , ...,x n ) 



MAS, 



(6) 



where 



d(S,u}i,u 2 , ■ ■ ■ ,Un-l) 
d(x 1 ,x 2 , ...,x n ) 



det 





df I 


dx\ 


dx n 
duji 


dx\ 


dx n 


duJn-l 


du n -i 


dx\ 


dx n _ 



(7) 



and oj\, . . . ,u; n -i are n — 1 solutions of Q or, equivalently, first integrals of ([5]) independent 
of each other. This means that M is a function of the variables (x\, . . . ,x n ) and depends on 
the chosen n — 1 solutions, in the sense that it varies as they vary. The essential properties of 
the Jacobi Last Multiplier are: 

(a) If one selects a different set of n — 1 independent solutions r/i,... ,rj n _i of equation fl3J), 
then the corresponding last multiplier N is linked to M by the relationship: 



N = M 



0(771 



,f/n-l) 



d(ui,. . .,u} n -i) 
(b) Given a non-singular transformation of variables 

t : (x±,X2, • • • , x n ) ► (x 1; x 2 , • • • , x n ), 
then the last multiplier M' of A'F = is given by: 



M' = M 



d(x\,X2, ...,X n ) 

d(x' l , x 2 , ■ ■ ■ , x' n ) 



where M obviously comes from the n—\ solutions of AF = which correspond to those 
chosen for A'F = through the inverse transformation t _1 . 
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(c) One can prove that each multiplier M is a solution of the following linear partial differ 
ential equation: 

d(M ai 



i=i 



dxi 



0; 



(8) 



viceversa every solution M of this equation is a Jacobi Last Multiplier. 



(d) If one knows two Jacobi Last Multipliers Mi and M2 of equation (Jl|) , then their ratio is a 
solution uj of dU), or, equivalently, a first integral of ([5]). Naturally the ratio may be quite 
trivial, namely a constant. Viceversa the product of a multiplier M\ times any solution 
uj yields another last multiplier M% = M\oj. 

Since the existence of a solution/first integral is consequent upon the existence of symmetry, an 
alternate formulation in terms of symmetries was provided by Lie |20j , [21] . A clear treatment 
of the formulation in terms of solutions/first integrals and symmetries is given by Bianchi [lj. 
If we know n — 1 symmetries of @/ ([5]), say 



,x n ) d x , i — 1 , h 1 , 



(9) 



Jacobi's last multiplier is given by M = A , provided that where 



det 



£1,1 



1.1 



Cl,n 
Cn— l,n 



(10) 



There is an obvious corollary to the results of Jacobi mentioned above. In the case that there 
exists a constant multiplier, the determinant is a first integral. This result is potentially very 
useful in the search for first integrals of systems of ordinary differential equations. In particular, 
if each component of the vector field of the equation of motion is missing the variable associated 
with that component, i.e., dai/dxi = 0, the last multiplier is a constant, and any other Jacobi 
Last Multiplier is a first integral. 

Another property of the Jacobi Last Multiplier is its (almost forgotten) relationship with 
the Lagrangian, L = L(t,u,u'), for any second-order equation 



u " = F(t,u,u') 



namely [30J: 



M 



d 2 L 
du' 2 



(11) 



(12) 
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where M = M(t, u, v!) satisfies the following equation 



dM dF 

■dT + M ft?= a (13) 

Then equation (fTT|) becomes the Euler-Lagrange equation: 

d (dL\ dL n 

The proof is given by taking the derivative of (|14p by v! and showing that this yields (|13|) . If 
one knows a Jacobi last multiplier, then L can be obtained by a double integration, i.e.: 

L = j (Jm&J\ &u' + f x {t,u)u' + f 2 {t,u), (15) 

where f\ and j% are functions of t and u which have to satisfy a single partial differential 
equation related to (fTTI) [27\. As it was shown in [27], /i, /2 are related to the gauge function 
g = g(t, u, u'). In fact, we may assume 

, dg 
ou 

h = % + h(t,u) (16) 

where has to satisfy the mentioned partial differential equation and g is obviously arbitrary. 
The importance of the gauge function should be stressed. In order to apply Noether's theorem 
correctly, one should not assume g = const, otherwise some first integrals may not be found 
(see [27] and the second-order equation in the next section). 

We now consider a fourth-order equation JT]) . In this case the Jacobi last multiplier satisfies 
the following equation 

dM dF 

■dT + M ^ = a (17) 

It is trivial to show from Fels' conditions j2])-([3]) that if a Lagrangian L = L(t, u, u',u") is taken 
such that 

M = ^ (18) 



along with the constraints^! 



du" 2 
dF 



du"' 

dF d / dF 



du 1 dt V du" 



0, (19) 
= (20) 



3 We note that (|20|) resembles the Euler-Lagrange equation (|14[) with F,u',u" replacing L,u,u' respectively. 
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then equation (pQ) becomes the Euler-Lagrange equation: 



^1 


(dL\ 




fdL\ 


dt 2 1 


\ du" I 


dt 1 


\du>) 



We emphasize that because of the assumption (|19p . then a Jacobi last multiplier M is trivial 
to find from equation (I17p . namely 

M = const. (22) 
Then L can be obtained by a double integration, i.e.: 

L = ^-u" 2 + h(t, u, u')u" + f 2 (t, u, u'), (23) 

where fi and f 2 are functions of t, u, v! which have to satisfy some partial differential equations 
related to ([T]). We can relate /i, f 2 to the gauge function g = g(t, u, u'). In fact, we may assume 



h 



dg_ 

du' 



h = ^u' + ^ + h(t,u,u') (24) 

where has to satisfy the mentioned partial differential equations and g is obviously arbitrary. 
Again we stress the importance of the gauge function. In order to apply Noether's theorem 
correctly, one should not assume g = const, otherwise some first integrals may not be found 
as in the examples of the next section. 

3 Two examples from Number Theory 

In this section we show in details the application of the Jacobi Last multiplier along with its 
connection to Lie symmetries. We also solve the inverse problem of calculus of variations as 
given in [6]. 

3.1 A second-order equation 

In [6] the following functional was introduced: 

y' A + 6uy 2 y' 2 dx, (25) 



o 



where y = y(x) £ C 2 [0,ir], y(0) = y(ir) = and v > 0. The corresponding Euler-Lagrange 
equation is 

y' 2 y" + vy 2 y" + vyy' 2 = . (26) 
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If we apply Lie group analysis to this equation, we find that it admits a two-dimensional abelian 
transitive Lie symmetry algebra (Type I) generated by the following two operators: 



yd y . 



(27) 



Then we can integrate equation (|26p . First, we introduce a basis of differential invariants of 
Ti, i.e.: 



u 



y_ 

y 



V = X . 



Then equation (|26p reduces to the following first-order equation: 

du —vuv 



dv uv 2 + u 2 ' 

which admits the operator T2 in the space of variables u, v, i.e. 

T 2 = vd v + ud u . 
Then its general solution is implicitly given by: 

u (2uv +u ) = const 



and in the original variable^: 



y' [2vy 2 + y /2 ) 1//4 = const 



viz 



vy 2 ai + ^Jv 2 y A a{ + 1) y 2 a\ 



yai \ uy 2 ai + y/ v 2 y* t a\ + 1 
with a\ an arbitrary constant Finally the general solution of (|26l) is given implicitly by: 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 



yai [vy 2 ai + ^Jv 2 y^a\ + 1 



vy 2 ai + ^v 2 y A d{ + l) y 2 ai 



d y = x + 02 



(34) 



We note that if y is positive and a\ = 1 the integral on the left-hand side could be integrated 
in terms of an hyper geometric function TC, namely: 



1 



2vy z H 





1 


1 






"1 


1" 


( 


~2' 


4' 


4_ 




2' 


2 



1 



y^v 2 



(35) 



4 The same first integral can be obtained by using Noether's theorem (see below) 
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Let us try to find a Lagrangian for equation f|26|) by using the Jacobi last multiplier, namely 
through (|12j) . The two Lie point symmetries (|27p yield a Jacobi last multiplier. In fact the 
following matrix [20], [2T] : 



/ i y' 
1 o 

\ y 



-y +y 



V 





2/ 



\ 



/ 



(36) 



has determinant different from zero and its inverse is a Jacobi last multiplier, i.e. 

2uy 2 + y' 2 



Mi 



y /2^ y 2 + y /2y 



(37) 



The corresponding Lagrangian is 



L 



1 

4i/y 



2i/y' arctan 
where f± , /2 are solutions of 



2vy 



+ log(2iV + y' 2 )i/y + 2 log(y>y ) + /i(x, y)y' + f 2 (x, y), 

(38) 

(39) 



0. 



dfi_dh 
dx dy 

If we impose the link between f%, f 2 with the gauge function g(x, y), namely (fT6|) . then f^(x, y) 
becomes just fz(x), an arbitrary function of the independent variable x. The Lagrangian (|38[) 
may appear ugly. Nevertheless the corresponding variational problem admits two Noether's 
symmetries, namely both Lie symmetries given in (|27p are Noether's symmetries. Consequently 
the following two first integrals of equation (|26p can be found by applying Noether's theorerd_|: 



Li 

r 2 



/i = (2iV + W, 



1 



dx + a>2 



Avy' 



2vy arctan 



2vy 



+ 2vy — Avxy 



(40) 

[g = s(x)y + x]. (41) 



with s(x) an arbitrary function of x. We note that the first integral I\ in (|40p was already 
derived in ([32]) . 

At this point one would like to know if it is possible to obtain the original Lagrangian given 
in ([25]), i.e. 

L H = y' 4 + 6vy 2 y' 2 (42) 

A property of the Jacobi last multiplier is that if one knows a Jacobi last multiplier and a first 
integral then their product gives another multiplier 26J. If we take the product of the first 



5 Also the corresponding gauge function g is given. It is important to remark that in the case of the first 
integral (|4ip the gauge function g cannot be constant, while it can be constant in the case of the first integral 
(1401) . Naturally, we have left out any inessential additive constants. 
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integral I\ P0|) and the multiplier M\ (|37|) . then we obtain another Jacobi last multiplier of 
equation ([26]) . i.e.: 

M 2 = -y' 2 - vy 2 (43) 
which can be integrated twice with respect to y' in order to yield the following Lagrangian^: 

L 2 = ~(y' 4 + Svy 2 y' 2 ) + hy 1 + h, (44) 

where f\ , f 2 are solutions of (|39|) . It is interesting to emphasize that this Lagrangian (namely 
Hall's Lagrangian) is such that the Lie operator T 2 in (|27p does not generate a Noether's 
symmetry for the corresponding variational problem. In fact only Ti is the generator of a 
Noether's symmetry for Hall's Lagrangian. 

If instead of (|25p we consider the functional [6j 

/ y' 4 + Qvy 2 y' 2 - 3A(i/)y 4 d x. (45) 
Jo 

and apply Lie group analysis to its corresponding Euler-Lagrange equation, viz: 

2/ 2/ + vy y +vyy + Xy = , (46) 

then we obtain the same Lie symmetry algebra generated by (|27|) which means that (|46p can 
be integrated by quadrature. In fact if we introduce the same variables as in (|28p then equation 
(|46|) reduces to the following first-order equation: 

du = -v {Xv 2 + vu 2 ) ^ 
dv u (yv 2 + u 2 ) 

which can be integrated to give 

Av 4 + 2vu 2 v 2 + n 4 = const (48) 

and in the original variables 

Ay 4 + 2zVV + y' 4 = const (49) 

a first integral of equation (|46p . 



Also Lie group analysis implies that if A = v 2 then equation (|46p admits an eight-dimensional 
Lie symmetry algebra generated by the following operators: 

Ti = —y (- cos(y/vx)d x + y^sin(^x)d y ) 
^2 = -y (sin(^x)d x + yy/v cos(y/vx)dy) 



Note the inessential multiplicative constant. 
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T3 = cos{2y/vx)d x — yy/v sm.{2y/vx)d y 

T4 = — sm{2\/ux)d x — y\/v cos{2\/vx)d y 

r 5 = d x (50) 

Te = yd y 

T7 = COs(\/vx)d y 

T 8 = — sm(y/vx)d y . 

which means that equation (I46|) is linearizable or indeed linear. Indeed in this case equation 
(JMD is just 

y" = -vy- (51) 

In order to find Lagrangians and first integrals for equation (j45]) we have to repeat mutatis 
mutandis the same analysis given above. We may underline that in the case of A = z/ 2 , namely 
equation (I51h . we can generate 14 different Lagrangians as it was shown in |27j . 

3.2 A fourth-order equation 

Another functional in [6] is the following: 

y' A + m 2 y" 2 dx. (52) 

The corresponding Euler-Lagrange equation is 

WlT + m 2 y m + 2 m ' 2 y" + 3^yy ,f2 - Qy' 2 y" = . (53) 

If we apply Lie group analysis to this equation, we find that it admits a three-dimensional Lie 
symmetry algebra generated by the following three operators: 

Xi = d x , X 2 = ydy, X 3 = xd x , (54) 

which means that we can reduce equation (|53p to a first-order equation, i.e.: 

du —7fiux — fiu — 6^x 3 — 4/ix 2 + Qx 
dx fiu 

with 







(55) 



- _ y"'y 2 2 y" 2 y 2 | y"y ~ _ y"y ^ 

If /x = 3 then equation f|53f) admits an eight-dimensional Lie symmetry algebra C generated by 
the following eight operators: 

3 x 3 
Ai = x 2 d x + -xydy, A 2 = xd x , A 3 = d x , A 4 = yd y , A 5 = — d y , 

^ y 

A 6 = —d y , A 7 = -d y , As = -d y (57) 

y y y 
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This means that equation (|53j) . i.e.: 

4y'y'" + yy m + 3y" 2 = (58) 

is linearizable by means of a point transformation [21 J . In order to find the linearizable transfor- 
mation we have to find an abelian intransitive two-dimensional subalgebra of C and, following 
Lie's classification of two-dimensional algebras in the real plane [21], we have to transform it 
into the canonical form 

d u , td u (59) 

with u and t the new dependent and independent variables, respectively. We found that one 
such subalgebra is that generated by A7 and As. Then we have to solve the following four 
linear partial differential equations of first order: 

A 7 (t) = 0, A 8 (i) = 0, A 7 {u) = t, A 8 (u) = l. (60) 

It is readily shown that the linearizable transformation is 

t = x, u = y 2 , (61) 

and equation (f58|) becomes 



u lv = 0. (62) 

Finally, the general solution of (|58|) is 

y = vVl + C12X + C13X 2 + a^x 3 (63) 
with a,i(i = 1,4) arbitrary constants. 

We note that if we apply the transformation (I6ip to equation (|53j) in the case of any fi then 
the following equation is obtained: 

iv _ (fi-3)(2uu" -u' 2 )u' 2 

which does not contain u'" , namely the third derivative of u by x. Therefore a constant, say 
1, is a Jacobi last multiplier of equation (f64"|) . and we can obtain a Lagrangian from ([23]) . i.e.: 

Li = ^u" 2 + h{x,u,u')u" + h(x,u,u') (65) 



where f± and satisfy (|24p and has to have the following expression: 

f 3 (x,u,u) = — — — 3 h hx{x,u)u' + h 2 (x,u) (66) 

24:fIU Z 
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with h\, hi arbitrary functions of x, u. 

We remark that a Jacobi last multiplier of equation (|53p can be derived from equation (|17l) . 

^-4^=0 -+ M-f, (67) 
M ax y 

although this Jacobi last multiplier is useless in order to find a Lagrangian of equation (|53p . 
Instead, if we apply transformation (|61p to (|65p in order to go back to the original function 
y(x), then adding also the particular assumptions fi = —u' 2 /(2u), and fa = u' 4 ([i + 1)/(16m 2 ) 
yield the following Lagrangian for equation (j53[) : 

Lmi = (m - I)?/ 4 + 2 2 / 2 / 2 (68) 



which, apart from an inessential multiplicative constant, is Hall's Lagrangian in (|52p . 

Let us apply Noether's theorem. If we consider the Lagrangian in (|65p . we find that the 

following two first integrals of equation (|64p can be obtained: 



+2imu'' i — fiu'^x — Quu' a + 3ii ,Zi x 
1 



X 3 + Xi =^ Ji = (12/MtV" - 4/iu 2 uy - 8fJ,u 2 u'u'"x + 4/iuV' 2 x 

'3 „,./4„ «„,„,/3 , 



X 2 =>■ J 2 = — ^(-fywtVu'" + 4/zuV 2 - nu' 4 + 3u' 4 ). (69) 
8[iu z 

A similar result is obtained if one uses Hall's Lagrangian and equation (|53p . Moreover, if we 
apply Noether's theorem to the linearizable equation (|58p with Hall's Lagrangian, i.e.: 

L H = y" + 3y 2 y" 2 , (70) 

we obtain the following seven first integrals^: 



Ai = 


> lmi 


= -3y 3 y" + 3y 3 y"'x + y 2 y' 2 + 7y 2 y'y"x 






-2y 2 y'y'"x 2 + y 2 y" 2 x 2 - 2yy' 3 x - Ayy l2 y"x 2 


-A 4 + A 2 = 


> Ini2 


= 9y 3 y"' + 2\y 2 y'y" - \2y 2 y'y"'x + 6y 2 y" 2 x 


A 3 = 


> I1TI3 


= -2y 2 y'y"' + y 2 y" 2 - Ayy' 2 y" + y" 


A 5 = 


> Irri5 


2 i n 1 // 2 1 III Z o/2 2 

= — 6y + byy x — iyy x + yy x — 6y x 


A 6 = 


> IrriQ 


= 2yy — 2yy" x + yy" x 2 — 2y' 2 x + 3y'y"x 2 


A 7 = 


> IrriY 


II , /// /2 , / // 

= -yy +yy x-y + 3yyx 


A 8 = 


> Im% 


in 1 / a 
= yy + 3y y 



/4„2 



f - 2Ayy' 2 y"x + 6y /4 x 



(71) 



7 Naturally, they are not all independent from each other. 
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Although we do not write down the corresponding expressions of the gauge function, we un- 
derline that it cannot always be set equal to a constant otherwise none of Imi, Im^, ImQ^mj, 
and img could be obtained. 

All seven first integrals (and even more) may be obtained without Noether's Theorem. 
In fact we just need to find the Jacobi last multipliers of equation (|58|) that are obtained 
by inverting the nonzero determinants of the possible 70 matrices made out of the eight Lie 
symmetries (I57D . Then the ratio of any two multipliers is a first integral of equation (j58|) . For 
example: 

W' + 3y" 2 \ 



C 



1234 



/ 


1 


y' 


y" 


y'" 




x 2 


\yx 


h - \ x y' 


y' ~ 




X 





-y' 


-2y 




1 











V 





y 


y' 


y" 



y 

-\y" - \xy"' 
-V" 



y'" 



(72) 



/ 



is the matrix obtained by considering the symmetries generated by operators Ai , A2 , A3 , A4 in 
(1571); its determinant is 



A 1234 = 9y"yy'y"' - iy%"' - \y l2 y" 2 + 6y" 3 y + \y"' 2 y\ 



and the corresponding Jacobi last multiplier is: 

1 



M 



1 



1234 



Ai2 3 4 I8y"yy'y>" - Ay' 3 y"' - 3y' 2 y" 2 + I2y" 3 y + 9y"' 2 y 2 ' 



(73) 



(74) 



Similarly, the matrix C5678 yields the determinant A5678 = 12/ y 4 , i.e. the Jacobi last multi- 
plieiH: 

M 56 78 = 2/ 4 (75) 

which we have already found in (|67p as an obvious solution of (|17p . Note that 



AT5678 



y\\%y"yy'y"' - ±y' 3 y"> - 2>y'^ + 12y" a y + 9y""y 



/2„//2 



,."3„ 



/"2„,2\ 



Mi2 34 

is "another" first integral of equation (1580 



(76) 



4 More fourth-order equations 

A similarity reduction of the fifth-order Korteweg-De Vries equation is the following fourth- 
order equation [2], [22] 



u 



20uu" + Wu' 2 - 40n 3 + au + t. 



(77) 



The multiplicative constant is inessential. 
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This equation satisfies Fels' conditions ©-([3]) as well as (fl~9]) - ([20]) . Therefore M = 1 is a Jacobi 
last multiplier, and formula (|18|) yields the following Lagrangian of (|77p : 

L = i(n" 2 + t 2 v!) + atuu - A0tu 3 u' + lOuu' 2 + ^c/(i, u, u') (78) 

where g = g(t, u, u') is the gauge function. 

In a paper on wave caustics [H], Kitaev derived the following equation: 



u iv = -^(Wuu'u'" -12au 5 + 8/3u 3 - 20u 7 + 40u V + 8u 4 t + 20u 3 u' 2 - 8uu' 
8u 2 

+I2uu" 2 - 8uu"t + uy 2 - \Qu' 2 u" + 8u' 2 t) (79) 

which satisfies Fels' conditions ©-(IS]), and admits M = u~ 2 as a Jacobi last multiplier. The 
simple transformation u = w 2 transforms equation (|79p into the following equation in w = w(t) 

w iv = —^(-12aw 8 + 8(3w A - 20iw 1 2 + 80wV + 160u> V 2 + 8w 6 t + 160WV 
low* 

-16W - I6ww"t - 80w' 4 + 16w' 2 t + t 2 ) (80) 

This equation satisfies Fels' conditions ©-([3]) as well as (fl~9]) - ([20j) . Therefore M = 1 is a Jacobi 
last multiplier, and formula (|18p yields the following Lagrangian for equation (|80p : 

L = - w" 2 + (-36<W^ + 24/3u/4 - 60w 12 t + 120u> V + 12w 6 t 2 + 40 W 3 
2 48w 3 

-2Aww't + t 3 ) + —g(t,w,w') (81) 
dt 

where 5 = g(t,w,w') is the gauge function. Replacing w = \pa into this Lagrangian yields the 
Lagrangian for equation (179]) . i.e.: 

L = — !-^-(12wV /2 - 12uu'V - 36au 5 u't + 24f3u 3 u't - 60uVt + 60uV 2 + 12uVi 2 
96tr 

-12W 2 t + nn't 3 + 8u 4 ) + -^(t, tt, u') (82) 

where 5 = g(t, u, u') is the gauge function. 

The method we proposed here is not the ultimate method. We can list many equations for 
which it does not work. In [23] seven fourth-order equations were derived as similarity reduc- 
tions of a mathematical model for thin liquid films [29] and its corresponding heir equations 
J . None of the seven equations passes condition ([3]) , although they all satisfy condition © . 
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5 Final remarks 



When one deals with a second-order differential equation the following remarks should be 
enlightened and kept in mind: 

• The most efficient Lagrangian, namely that which allows the most number of Noether's 
symmetries, may not be the Lagrangian with the simplest form. 

• Lie symmetries are the key tool for finding Jacobi last multipliers and therefore La- 
grangians. 

In [5] the necessary and sufficient conditions under which a fourth-order equation ([1]) admits 
a unique Lagrangian were determined. In this paper we suggest that the application of the 
Jacobi last multiplier may yield that unique Lagrangian in some instances. 

If one rewrites a fourth-order equation as either a suitable system of two second-order 
equations [4], or a system of four first-order equations [8] the challenge of solving an inverse 
problem of calculus of variations may still be open. 
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